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Summary. Self-regulation of living tissue as an example of self-organization phe- 
nomena in hierarchical systems of biological, ecological, and social nature is under 
consideration. The characteristic feature of these systems is the absence of any gov- 
erning center and, thereby, their self-regulation is based on a cooperative interaction 
of all the elements. The work develops a mathematical theory of a vascular network 
response to local effects on scales of individual units of peripheral circulation. 



1 Introduction: Living Tissue Self-Regulation Problem 

Biochemical compounds (oxygen, nutrients, etc.) required by cells of human 
body for their life are transported by blood flow through the vascular network. 
Simultaneously, blood flow withdraws products of their life activity, e.g., car- 
bon dioxide. So when cells in a certain region small in comparison with the 
whole organism consume an increased amount of oxygen and nutrients, blood 
flow through it just has to grow to maintain the internal environment of the 
cellular tissue. This effect will be referred to as the vascular network response 
to local variations in the state of the cellular tissue. 

We consider the response of the vascular network to be perfect if every 
cell is supplied with such amount of oxygen and nutrients that is strictly 
necessary for its current activity. As a result, if the functioning of a certain 
tissue region has changed and cells in it need an increased amount of oxygen 
and nutrients, then the vascular network with the perfect behavior is able to 
deliver the required biochemical compounds to this region in such a manner 
that blood flow at the other points of the cellular tissue not be disturbed 
at all. However, in the general case, the vascular network responding to the 
local need of cells at one point will disturb the functioning of cellular tissue 
at all its other points. This effect can give rise to an erratic behavior of the 
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organ as a whole and, thus, should be depressed. Since there are no special 
centers controlling blood flow redistribution on scales of individual organs this 
regulation can be implemented only via a cooperative interaction between all 
its elements. The latter feature allows us to call such regulation processes 
distributed self-regulation. 

Up to now specific mechanisms by which the distributed self-regulation 
arises are far from being understood well. In order to construct an appropriate 
model for the distributed self-regulation two aspects must be elucidated. The 
first one is the self-processing of information because none of the blood ves- 
sels possesses the complete information about the state of the corresponding 
organ. Moreover, each vessel really needs only a piece of information required 
for its individual functioning. The second aspect concerns the cooperative 
mechanism governing the blood flow redistribution over the vascular network. 
Indeed, again, none of the vessels controls individually blood flow even through 
itself; responding to the appropriate information it can only dilate or constrict 
to change its hydrodynamic resistance to blood flow. 

Here we present a model for the semi-perfect self-regulation of living tis- 
sue that does not require the vessels response to be ideal as was presumed 
previously [1, 2], which enables us to describe observed phenomena in the real 
living tissue. The key point of the proposed model is the reaction of the blood 
vessels to the concentration of special biochemical components, activators, 
generated by the cells during their life activity most intensively when one of 
the internal environment parameters comes close to the tolerance boundary. 
The concept of activators can be justified appealing to biomedical data [3]. 

2 Living Tissue Model 

The model for the vascular network architecture assumes the following. First, 
the arterial and venous beds are individually of the tree form and the branch- 
ing of vessels is symmetrical. Second, the vessels of each level are distributed 
uniformly in the living tissue domain considered to be a square or cube. The 
first assumption enables us to introduce the vessel hierarchy and to order all 
the arteries and veins according to their position in the hierarchical structure. 
The second one allows us to specify the vessel arrangement using a self-similar 
embedding illustrated in Fig. 1 . Besides, for the sake of simplicity, the arterial 
and venous beds are assumed to be the mirror image of each other with respect 
to all the properties, including the vessel response. This artificial assumption 
just simplifies the description of the vascular network response. 

Blood flow distribution over the vein network {v} (Fig. 1) is described 
as follows. First, to every vein v we ascribe the blood flow rate J v and the 
activator concentration V in blood going through it. The state of this vein is 
quantified in terms of the vessel resistance R n {0 v ) depending on the activator 
concentration V , which reflects the vessel response. Here the function R n (9) 
is assumed to be identical for all the veins of one level n. Second, to every 
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Fig. 1. Vascular network model reduced to the venous bed and its reaction to 
the distribution of activators in the cellular tissue. Here, in particular, Q — Q c 
designates the elementary domain related to the venule c. Left fragments illustrate 
the vessel embedding into the cellular tissue. 



branching point k we ascribe a certain blood pressure P. Then the blood flow 
distribution is governed by the system of equations including Poiseuille's law 
written for every vein v 

Pin — -Pout — JvRn{&v) , (1) 

and conservation of blood and activator substance at every branching point k 

JinA + J'm,2 = Jout , (2) 
•^in,l^in,l + </in,2#in,2 = ^out^out ■ (3) 

At the entrance to the venules {c}, i.e. the veins of the last level iV and the 
exit from the root vein (the vein of level n = 0) 

P\n\c = P > Pout|n=0 = . (4) 

Finally, the equality of the activator concentration in every venule c and the 
mean concentration of activators in the corresponding elementary domain Q c 

0c = B{Qc) (5) 

holds. The elementary domains {Q} form the partition of the living tissue bulk 
in such way that every elementary domain contains just one venule, which is 



4 



Wassily Lubashevsky, Ihor Lubashevsky, and Reinhard Mahnke 



illustrated in Fig. 1. In the given analysis the distribution of activators in the 
cellular tissue specified in terms of {9(Q)} is treated as given beforehand. 

When one of the internal environment parameters comes close to the 
boundary of tolerance zone, the living tissue to survive has to increase es- 
sentially the blood flow rate via vessel dilation. So blood vessels should ex- 
haust their ability to widen when the activator concentration exceeds a certain 
threshold A matching these critical conditions. It allows us to adopt the fol- 
lowing ansatz 

Rn(0v) = Pn<t>(0v) , (6) 

where p n is the vessel resistance at 9 V = which depends only on the number 
of the hierarchy level and (f)(9) is a certain function universal for all the vessels. 
By definition, at 9 = the equality (f)(0) = 1 holds and for 9 » 4 the 
inequality (f)(9) — > (f>n m <C 1 should be the case. The following ansatz 
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is used to be specific. Here the additional numeric parameter e characterizes 
the linearity of the function (f)(9) in the region 9 < A. Figure 2 visualizes 
the given function (f)(9) for various values of its parameters. The values e = 
and </>i; m = match the ideal vessel response. Ansatz (7) can be justified, for 
example, appealing directly to Fig. 2. 
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Fig. 2. The vessel response function 
of its parameters. 



specified by ansatz (7) for various values 



The dependence of the vessel resistance on the hierarchy level, i.e. the 
value of p n vs the level number n is constructed using the known equality of 
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contributions to the vascular resistance of arterial beds from all the hierarchy 
levels. This enables us to write the /?„-dependence in the form 

Pn = (20 n po , (9) 

where the value po is related with the root vein and the parameter £ 1 allows 
us to consider the cases where various groups of vessels dominate in the blood 
pressure distribution. In particular, for £ < 1 the large veins contribute mainly 
to the vascular resistance of the venous bed, whereas, in the opposite case, 
C > 1, the small veins are dominating. 

3 Results 




Fig. 3. Example of the systems analyzed numerically. The region with a nonzero 
value of the activator concentration 9 is shadowed and the ordering of the venules 
is illustrated for the system with the number of the hierarchy levels equal to N = 6. 

The blood flow redistribution governed by equations (l)-(3) subjected to 
conditions (4) and (5) was studied numerically. Let us, first, present the results 
for the case when the activator concentration 9 differs from zero inside m first 
elementary domains located at one of the square corner (2D model of living 
tissue, Fig. 3), namely, 

Oc.i = if 1 < i < m and 9 Cji = if 1< m < N . (10) 

Figure 3 illustrates the analyzed case and shows the introduced ordering of 
the venules with index i. The elementary domains are also labeled with the 
index i of the corresponding venule. 

Figure 4 presents the obtained results. The first column exhibits the de- 
pendence of the blood flow rate J\ through the venule v\ on the activator 
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Fig. 4. Results of numerical analysis. The blood flow rate Ji (in units of In) vs 
the concentration of activators in the elementary domain Qi (left column) and the 
distribution of the blood flow rates (in units of In) over the venules (right column). 
In obtaining these results the parameters e — 1.0 and </>n m = 0.1 were used, the used 
values of the other parameters are shown in figure. 

concentration 9 in units of the blood flow rate In through venules under the 
normal conditions, i.e. in the case of absence of activators. The number of 
the hierarchy levels was set equal to N — 6; it meets tenfold decrease in the 
vessels length when passing from the root vessel to the smallest arterioles or 
venules. 

The left upper frame depicts three curves corresponding to different values 
of the parameter (. When the activator concentration in the domain Qi gets 
the threshold 9 C — A of the individual vessel reaction the venule joined with 
this domain exhausts its capacity to widen. However it is not enough for the 
blood flow rate to grow substantially. Indeed, the contribution of this venule 
to the vascular resistance is small even along the path Ve leading from the 
root vessel to the domain Qi on the vein tree (Fig. 3). So many veins along 
this path have to dilate in order for the blood flow rate to grow essentially in 
Qi. However, the larger a vein, the larger the tissue region drained through 
this vein as a whole and its size increases exponentially with the number of 
branching points separating the venule v\ and the given vein along Vg. As 
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a result when the activator concentration in the domain Q x is about 9 ~ A 
the activator concentration in relatively large veins on the path Vg turns 
out to be rather small and these vessels cannot widen remarkably. Thereby 
the blood flow rate cannot also exhibit substantial increase (Fig. 4). As the 
activator concentration grows further the induced increase in the activator 
concentration along the path Vg gets the critical value A in larger and larger 
veins, causing the growth of the blood flow rate Ji . When the main part of 
the veins belonging to Vg widen to the maximum the blood flow rate comes 
to the upper limit J max . 

The maximum J max of the blood flow rate depends on the system param- 
eter, in particular, on £. As it will be demonstrated below practically all the 
additional amount of blood entering the system due to the vessel dilation is 
directed to the domain Q\. So the higher the contribution of large veins, the 
higher the upper limit of the blood flow rate. These speculations are justified 
well in Fig. 4, the left upper frame. Because of this blood flow focusing along 
the path Vg the value J max can exceed substantially the value l/0n m which 
would be attained if the total blood pressure drop had fallen just on the last 
vessel supplying the region Qi with blood. A similar effect should be expected 
when the region of living tissue where activators are located becomes larger. 
So the value of J ma x has to drop with m increasing, which is demonstrated in 
Fig. 4 (left bottom frame) . The value J max goes down to 0n m when activators 
spread over the whole region of the microcirculatory bed. 

The obtained results demonstrate us also the fact that when the vessel 
behavior is not ideal the threshold 6 t h in the vascular network response as 
a whole does not coincide with the critical value A of the individual vessel 
reaction and can exceed it remarkably. Let us treat the threshold 9th as the 
value of 9 at which the blood flow rate attains one half of its upper limit. 
Then this fact is clearly visible in Fig. 4 (left upper frame) and, for example, 
the estimate 0th ~ 2Z\ holds for the case shown by curve 1. Naturally, the 
threshold 0th depends also on the parameters of the vascular network but not 
only on the characteristics of individual vessel behavior. 

The shown curves enable us to assume that due to the hierarchical struc- 
ture of the vascular network there is a certain function r){6\<j>, e) := J /In 
depending only on the parameters of individual vessel behavior, which de- 
scribes the increase in the blood flow rate with the growth of the activator 
concentration until the latter becomes greater then the threshold t h- The 
other characteristics of the vascular network and the distribution of activa- 
tors over the cellular tissue, roughly speaking, affect only the upper limit J max 
of blood flow rate. 

The right column depicts the distribution of the blood flow rates over 
the venules. Comparing the data shown hear and that of the left column we 
can declare that the main increase in the blood flow rate is located along 
the path Vg. However, the question about the interference of blood streams 
induced by distant regions of living tissue requires an individual investigation. 
Nevertheless, when the living tissue is excited only in one domain its response 
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can be treated as quasilocal. We have used the term "quasilocal" because this 
response also depends on the size of the excited domain and contains a certain 
component describing small variations in the blood flow rate at distant points. 

As the nonlocal component of the Ji{6>}-functional is concerned, the right 
column of Fig. 4 demonstrates us that depending on the relation between 
the activator concentration 9 in the elementary domain Qi and the vascular 
network threshold 9th it changes its form. When 9 < 6th the increase in the 
blood flow rate through Qi is caused by the blood flow decrease through 
the other elementary domains. In this case the required growth of the blood 
flow rate is due to the redistributed of blood flow over the vascular network 
without additional amount of blood entering the microcirculatory bed. Since 
the number of the surrounding domains is rather large the decrease in the 
blood flow rate outside the excited region is not essential. When 9 > 6 t h the 
increase in the blood flow rate is mainly caused by additional portion of blood 
entering the microcirculatory bed. Therefore the blood flow rate is increased 
in the surrounding domains also. 




Fig. 5. Results of numerical analysis. The family of the dependencies {Ji{0i)} each 
of them is the rate Ji of blood flow through venule i vs the activator concentration 
6i in it for the activator distribution given by expression (11). In obtaining these 
results the parameters e = 1.0 and </>ii m = 0.1 were used, the used values of the other 
parameters are shown in figure. 



To analyze the validity of the quasi-locality approximation of the living 
tissue response we have studied the blood flow rate distribution Ji over the 
venules for the exponential form of the activator distribution over the cellular 
tissue, namely, 

(11) 



#exp 



(i-1) 



Living Tissue Self-Regulation Problem 



9 



for different values of the parameter a. Figure 5 exhibits the resulting depen- 
dence of the rate Ji(9i) of blood flow through venule i on the local value Oi 
of the activator concentration in this venule. Plotting the family of the curves 
Ji(9i) for a collection of venules we can verify whether the dependence Ji(9i) is 
really a local function of its argument. The found results (Fig. 5) demonstrate 
us that when the activators are concentrated in the region whose size is larger 
then the size of elementary domain the blood flow rate Ji does depend only 
on the local value of the activator concentration 9i. Only when the region of 
activator localization is comparable with the elementary domains in size the 
quasi-locality approximation is violated to some degree. 

4 Conclusion 

We have presented a model for the distributed self-regulation of living tissue 
that exemplifies a general mechanism by which self-regulation of active hier- 
archical systems without controlling centers can arise. The key point is the 
information self-processing caused by conservation of blood and activators in 
blood stream through the arterial and venous beds. It is essential that the 
vascular network response is based on vessels reacting to certain biochemical 
components, activators, that are generated by the cells when one of the in- 
ternal environment parameters comes close to the tolerance boundary. In this 
case under rather general assumptions about the specific properties of this 
reaction the vascular network response to local perturbations is quasi-local. 
In other words the quasi-locality is the consequence of the hierarchical struc- 
ture of the vascular network rather then the ideal behavior of its elements as 
though previously [1, 2]. In this way we have overcome the main obstacle to 
applying the developed theory of distributed self-regulation to describing real 
phenomena in such systems. Besides, due to the hierarchical structure of the 
vascular network its response turns out to be much stronger then it could be 
expected considering the individual dilation of blood vessels. 
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